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abstract 



In this paper, we investigate the evolution of classical wave propagation in the canonical acoustic 
black hole by numerical method and discuss the details of tail phenomenon. The oscillating frequency 
and damping time scale both increase with the angular momentum For the lower I, numerical 
O I results show the lowest WKB approximation gives the most reliable result. We also find that time 

scale of the interim region from ringing to tail is not affected obviously by changing I. 



I. INTRODUCTION 



Some properties of black holes can be investigated using acoustic analogues in the laboratory through the propa- 
gation of sound wave. Hawking radiation is a remarkable prediction and is almost universally believed to be one of 
the most important in black hole's physics. However, the Hawking temperature of astrophysical black holes is much 
smaller than the temperature of the cosmic microwave background so that one cannot acquire any conclusive evidence 
of the existence of Hawking radiation. About twenty-five years ago, Unruh proposed a method that certain aspects 
of astrophysical black hole are mapped into problems in the theory of supersonic acoustic flows [l[ . Even though the 
Hawking temperatures associated to acoustic analogues are not high enough to be detectable up to now, the situation 
is likely to change in the near future .2]. A profound understanding of the classical physics of acoustic black hole is 
. indispensable for the detection of Hawking radiation. Berti, Cardoso and Lemos 0] investigated wave propagation in 
the "draining bathtub" model and the "canonical" (l+3)-dimensional acoustic black hole Especially, using the 
Wentzel-Kramers-Brillouin (WKB) method, they calculated the quasinormal modes (QNMs). Many physicists believe 
that the figure of QNMs is a significant fingerprint indirectly identifying the existence of a black hole. The QNMs of 
black holes in the framework of general relativity [j| @ and string theory @, § has been studied widely. 

Approximately, there are three stages in the evolution of the perturbations of an acoustic black hole [j] . First stage 
is the rapid response at very early time, on which the initial conditions have a great effect. Second stage is quasinormal 
ringing phase, which characteristic oscillation frequencies and damping times depend strongly on the acoustic analogue 
QNMs. The QNMs are determined completely by the parameters of system, therefore they would carry significant 
information about the background curvature of the intervening spacetime. Finally, there is a tail stage, which decays 
approximately as a power in time owing to backscattering off the spacetime curvature. In Ref. 3, the authors have 
used three WKB computational schemes, i.e. the lowest approximation 0, 3rd order improvements [13, El and 6th 
order corrections For the canonical acoustic black hole, the results show that / = 1 QN frequencies seem to 

be the problem, in which the mode suffers a large variation as one goes from the lower approximation to the higher 
approximation. This means that the WKB approach is more dependable for higher I, which was first discovered in 
the early work [1, [HJ [Hj] . Therefore, the numerical calculation is necessary for the lower I QN frequencies. 

In this paper, we investigate in detail the relations between QNMs of canonical acoustic black hole and the angular 
momentum I by the numerical calculation in null coordinates. Some results attained by this way are supported by 
the analytic results and WKB results. Most of importance, we confirm that the lowest WKB approximation gives 
the most reliable results for I = 1 case. Furthermore, we show a picture of classical wave propagation including the 
interim region from the quasinormal ringing to tail stage. 
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II. FORMALISM AND BASIC EQUATIONS 

Assume the fluid to be incompressible and spherically symmetric, then since background density p is position 
independent the continuity equation implies the velocity is in proportion to r~ 2 . The background pressure p and 

speed of sound c are also position independent because of the barotropic assumption. Therefore, one can define a 

i 

normalization constant r*o = (vr 2 /cj 2 . The canonical acoustic metric describing the propagation of sound waves in 
this incompressible and spherically symmetric (l+3)-dimensional fluid flow [J] is: 

ds 2 = -c 2 - ^ dt 2 + (l - $\ 1 dr 2 + r 2 (d9 2 + sin 2 0d<f> 2 ) (1) 

The metric ([I} is distinct from any of the geometries typically considered in general relativity. Unruh [l[ first suggested 
that the propagation of a sound wave is described by the Klein-Gordon equation V^V^"]/ = for a massless scalar 
field ^ in a Lorentzian acoustic geometry, which take metric QJ at present. We can separate variables by setting 

*(t,r,M) = -*(r.,*)3WM) (2) 
r 

where Yi m {9, <f>) are the usual spherical harmonics and the tortoise coordinate r* is defined by 

r*=/(l- (3) 
where we have chosen unit c = 1. The evolution equation of $(r*) is 

<9 2 $ <9 2 $ 

-^ + ^r v " (4) 

where the effective potential 

rw-(i-2)[!S£> + £] (5) 

We introduce the null coordinates u — t — and w = £ + r* , Eq. ([5]) can be reduced to 

<9 2 $ 

- 4 ^ = y( ^ (6) 
Eq. ([5]) can be numerically integrated by the ordinary finite element method. Using the Taylor expansion, we have 

§n = *b + <S>w - *s - SuSv(— — — — — -)^— — -V(u) + 0(A 4 ) (7) 

where N, W , £J and S 1 are the points of a unit grid on the u — v plane which correspond to (u + A,i> + A ),{u + 
A,w),(u,u + A) and (it, u), and A is the step length of the change of u or v, i.e., A = 5u = 5v @, 0. Because 
the quasinormal ringing stage and the late time stage are both insensitive to the initial conditions, we begin with a 
Gaussian pulse of width a centred on v c when u — u$ and set the field $ is zero on v — vq, 

®(u = u ,v) = exp[ 2 ^ 2 C ] 

$(«,w = «b) = (8) 

Next, the point in the u — v plane can be calculated by using Eq. @, successively. Finally, the values of <&{u max , v) 
are extracted after the integration is completed where u max represents the maximum of u. Taking sufficiently large 
Umax for the various w-value, we obtain a good approximation for the wavefunction of canonical acoustic black hole. 
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III. NUMERICAL RESULTS 



Our numerical results, which are all consistent with the analytic results and WKB results in Ref. 3, are shown in 
Figs. 1-6. As a reminder, the oscillating period, damping time scale and late time tail are shown in these figures. 
Here, the parameter r is set to uninty. The dependence of quasinormal modes on r is trivial. On the one hand, 
the canonical acoustic metric coordinates can be rescaled to set r = 1. On the other hand, the results must depend 
linearly on r since it is the only dimensional quantity in the problem. In Fig. 1, we show the relations between 
the wavefunction and the angular momentum I. Our numerical result is consistent with Ref. 3. That means the 
oscillating period and the damping frequency both decrease when the index I increases. Furthermore, we confirm that 
the lowest WKB approximation gives the most reliable results for I = 1 case. To further corroborate this conclusion, 
we list QN frequencies for I = 1,2,3,4 in Table 1. 



QN frequencies of the canonical acoustic black hole for I — 1, 2, 3, 4 are listed. 



1 


Re(tu) 


I in (u>) 


1 


1.463 


0.666 


2 


1.619 


0.653 


3 


1.642 


0.625 


4 


1.758 


0.620 



In Figs. 2-6, we choose I = 2, and consider in detail the picture of classical wave propagation in the canonical 
three stages, the second and final stages as illustrated in Fig. 2. The prompt contribution is the evident counterpart 
of light cone propagation in the V = case, which strongly depends on the initial conditions, therefore it is left 
out in Fig. 2. At intermediate v values the wavefunction is dominated by an exponential decay, whose oscillation 
frequency and damping time are described by its QNMs. At the late-time (large v value) the propagating wave leaves 
a power-law tail which is magnified in Fig. 3. By numerical calculation, we attain the expression of power-law falloff, 
$ w 7.36 x 10~ 32 i~ 10 , which is consistent with the analytic result in Ref. 3. Especially, the interim region from 
ringing stage to tail stage is corresponding to rectangular region A in Fig. 2. This interim region is replotted in 
Fig. 4, where the w-coordinate is magnified about 10 3 times. The time interval from region B to region C is so short 
that the numerical results between region B and region C seem unfaithful. Therefore, we do not discuss physical 
implications about this region, attentively. The rectangular regions B and C of Fig. 4 are magnified in Fig. 5 and 
Fig. 6, respectively. In region B, the oscillation frequency dramatically changes and tends to zero. Likewise, in region 
C, the damping time scale also has a drastic change, which becomes infinity. These figures tell us how perturbation in 
vicinities of this black hole die out as a late-time tail. It is easy to find that time scale of interim region from ringing 
to tail is not affected obviously by changing the angular momentum I. 



IV. CONCLUSIONS 



In this work we considered numerically the evolution of classical wave propagation in the canonical acoustic black 
hole and discussed the details of tail phenomenon. We summarize main results as follows: 




FIG. 1: For the canonical acoustic black hole with ro = 1, the wavefunctions are shown via different the angular momentum I. 
The logarithm is to base 10. 
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(i) For I > 2, the numerical results are consistent with the first, third and sixth order WKB method. For the 
lower I, numerical results show the first WKB approximation gives the most reliable result because of the basic WKB 
assumption (the ratio of the derivatives of the potential to the potential itself should be small) is broken. 

(ii) From a physical viewpoint, the most reasonable explanation for the production of late-time tails is the backscat- 
tering of waves off a spacetime curvature at asymptotically far regions. Our numerical results show that late-time 
tail is consistent with Rcf. 13, and time scale of the interim region from ringing to tail is not affected obviously by 
changing the angular momentum I. 

(iii) The oscillating frequency and damping time scale both increase with the angular momentum /. In the limit of 
large I, the real part of fundamental QN frequency increase linearly and imaginary part tend to a constant with the 
angular momentum I. 
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FIG. 2: For the canonical acoustic black hole with r = 1 and I — 2, the rectangular region A describes as interim region from 
ringing to tail. The logarithm is to base 10. 
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FIG. 4: For the canonical acoustic black hole with ro = 1 and I = 2, the interim region from ringing to tail is shown. 




FIG. 5: The rectangular region B in Fig. 4 is magnified. The logarithm is to base 10. 




FIG. 6: The rectangular region C in Fig. 4 is magnified. The logarithm is to base 10. 




FIG. 3: The late-time tail 



in Fig. 2 is magnified. The logarithm is to base 10. 



